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ABSTRACT 


Mimulus guttatus DC. (yellow monkey-flower; Phrymaceae) is an important model species for ecological and evolutionary 
studies, being locally adapted to a wide range of elevation, moisture and temperature gradients, soil types, and pollinator 
availabilities. In order to advance this species as a model for evolutionary genetic studies, we have developed virus-induced 
gene silencing (VIGS) using the tobacco rattle virus (TRV) to assay gene function. We demonstrate the effectiveness of 
Agrobactertum-mediated VIGS in two divergent populations of M. guttatus, Iron Mountain 767 (IM767) and Point Reyes (PR). 
Plants infected with a fragment of the carotenoid biosynthesis pathway gene PHYTOENE DESATURASE (PDS) cloned into the 
TRV2 vector exhibited endogenous PDS silencing and photobleached phenotypes. We further assayed for VIGS-induced floral 
phenotypes by silencing paralogous genes putatively affecting floral symmetry, CYCLOIDEAI (CYC1) and CYCLOIDEA2 
(CYC2). Simultaneous silencing of CYC] and CYC2 resulted in organ number defects in the petal and stamen whorls; silencing 
of CYC] affected petal margin growth; and silencing of CYC2 had no effect on flower development. Infection with TRV2 and 
TRV1 is significantly higher and more pervasive in the IM767 versus the PR population and is more efficient after vacuum 
infiltration. These results demonstrate the efficacy of VIGS for determining the function of developmental genes, including those 


involved in ecologically important reproductive traits. 
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The cosmopolitan plant genus Mimulus L. (Phry- 
maceae) comprises ca. 160 to 200 species, the 
majority of which are derived from two rapid 
radiations in western North America and Australia 
(Beardsley & Olmstead, 2002). Species of Mimulus 
occupy a vast range of habitats, from desert to alpine, 
and are phenotypically diverse in both reproductive 
and life history traits (Vickery, 1978; Beardsley et al., 
2004; Angert & Schemske, 2005; Wu et al., 2008). 
For example, major differences in floral form are 
associated with multiple shifts in pollinators and 
changes in the relative importance of selfing versus 
outcrossing (Bradshaw & Schemske, 2003). Much of 
the ecological and morphological diversity within the 
genus is evident in the young, but genetically diverse, 
species complex M. guttatus DC. (yellow monkey- 
flower) (Kelly & Willis, 1998; Sweigart & Willis, 
2003; Beardsley et al., 2004). Populations of the M. 
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guttatus complex are locally adapted to a range of 
elevation, moisture and temperature gradients, soil 
types, and pollinator availability, making them an 
excellent system for ecological studies (Allen & 
Sheppard, 1971; Fishman et al., 2002; Hall & Willis, 
2005; Wu et al., 2008). Recent development of 
genomic tools, including a near-complete genome 
sequence of M. guttatus, have further made Mimulus 
an exciting model group in which to test fundamental 
evolutionary questions related to speciation, adapta- 
tion, and macroevolution (Hall & Willis, 2006; 
Holeski & Kelly, 2006; Cooley & Willis, 2009; Ivey 
et al., 2009; Lowry et al., 2009; Wu et al., 2010). 
In addition to genomic tools, the development of 
Mimulus guttatus for evolutionary developmental 
genetic (“evo-devo”) studies requires transformation- 
al capabilities to assay gene function. This can either 
be accomplished through forward or reverse genetic 
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approaches. One promising reverse genetic method 
for the rapid and transient (intragenerational) silenc- 
ing of gene expression is virus-induced gene 
silencing (VIGS) (Baulcombe, 1999; Robertson, 
2004; Purkayastha & Dasgupta, 2009; reviewed in 
Di Stilio, 2011). VIGS utilizes the RNA interference 
(RNAi) pathway of plants to target endogenous 
messenger RNA (mRNA) for degradation (Fagard & 
Vaucheret, 2000; Robertson, 2004). Specifically, 
Agrobacterium-mediated infection of plant cells with 
a modified viral vector, containing a plant gene of 
interest, results in the expression of double-stranded 
RNAs that are degraded in vivo by DICER-like 
enzymes (Hamilton & Baulcombe, 1999; Zamore et 
al., 2000). This results in short, interfering RNAs of 
12 to 23 bp that form the template for further 
degradation of complimentary RNAs, including 
endogenous mRNAs (Hammond et al., 2000; Zamore 
et al., 2000). 

VIGS has proven effective in a wide range of 
distantly related angiosperm taxa (Di Stilio, 2011), 
including the monocots Brachypodium distachyon 
(L.) P. Beauv. (purple false brome) (Demircan & 
Akkaya, 2010), Hordeum vulgare L. (barley) (Holz- 
berg et al., 2002), orchids (Lu et al., 2007), Triticum 
aestivum L. (wheat) (Scofield et al., 2005; Tai et al., 
2005; Scofield & Nelson, 2009), and Zingiber 
officinale Roscoe (culinary ginger) (Renner et al., 
2009); the basal eudicots Aquilegia vulgaris L. 
(common columbine) (Gould & Kramer, 2007), 
Eschscholzia californica Cham. (California poppy) 
(Wege et al., 2007), Papaver somniferum L. (opium 
poppy) (Hileman et al., 2005), and Thalictrum sp. L. 
(meadow rue) (Di Stilio et al., 2010); the asterids 
Antirrhinum majus L. (snapdragon) (Preston & Hile- 
man, 2010), Nicotiana benthamiana Domin (tobacco) 
(Kumagai et al., 1995), Petunia hybrida (Hook.) Vilm 
(petunia) (Spitzer et al., 2007), Solanum lycopersicum 
L. (tomato) (Liu et al., 2002), and Spinacia oleracea 
L. (spinach) (Sather et al., 2010); and the rosids 
Arabidopsis thaliana (L.) Heynh. (thale cress) (Burch- 
Smith et al., 2006; Pflieger et al., 2008; Igarashi et 
al., 2009), Glycine max (L.) Merr. (soybean) (Zhang et 
al., 2009), and Medicago truncatula Gaertn. 
(Grønlund et al., 2008). However, despite the 
generality of VIGS, it is clear that not all vectors 
are effective in all species and that protocol 
optimization is required to elicit efficient gene 
silencing in nonmodel taxa (Liu & Page, 2008; Di 
Stilio, 2011). 

Here, we demonstrate the effective use of VIGS in 
two morphologically diverse populations, Point Reyes 
(PR) and Iron Mountain 767 (IM767), of the emerging 


model species Mimulus guttatus using constructs 
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containing PHYTOENE DESATURASE (PDS) or 
paralogous CYCLOIDEA-like (CYC-like) genes. PDS 
is a functionally conserved plant gene that is involved 
in the carotenoid biosynthesis pathway that protects 
the chloroplast from photo-oxidation; PDS silencing 
causes a characteristic photobleaching phenotype 
under high light growth conditions (Demmig-Adams 
& Adams, 1992; Benedito et al., 2004; Hileman et 
al., 2005; Gould & Kramer, 2007). By contrast, CYC- 
like genes generally affect cell proliferation and 
expansion in the petal and stamen whorls, mainly in 
the dorsal region of the flower (Luo et al., 1996, 1999; 
Cubas et al., 1999; Feng et al., 2006; Busch & 
Zachgo, 2007; Broholm et al., 2008; Fambrini et al., 
2011). 


MATERIALS AND METHODS 
CONSTRUCTION OF MIMULUS GUTTATUS TRV2 VIGS VECTOR 


RNA was extracted from Mimulus guttatus PR and 
IM767 material using the RNeasy Plant Mini Kit and 
RNase-Free DNase Set (Qiagen, Valencia, Califor- 
nia), and complementary DNA (cDNA) was generated 
using the iScript cDNA Synthesis Kit (Bio-Rad, 
Hercules, California). MgPDS was amplified from 
cDNA of both M. guttatus populations using the 
conserved primers ALL.PDS-F (5’-TTGWATGC- 
CAARYAARCCWG-3’) and ALL.PDS-R (5’-ACAC- 
TIAGMAGKGRRCTTCTGC-3"), which carry degen- 
eracies in order to amplify PDS from a number of taxa 
in the Lamiales. MgCYC] and MgCYC2 were 
amplified from IM767 cDNA using the primers 
MeCYC1.full.F (5’-ACATGATTCCTCATGGCTCTT- 
CTC-3") and MgCYCl1.full.R (5’-ACTTGTGGT- 
GATCCAAAATGGCAG-3"), and MgCYC2.full.F 
(5'-CAGTTACCTCCATCCTCAGGCAAC-3") and 
MgCYC2.full.R (5’-ACTTGTGCTGGTCCAAAATA- 
GCAC-3’). The resulting fragments were cloned 
separately into TOPO-TA 2.1 (Invitrogen, Carlsbad, 
California), and their identity was confirmed by 
sequencing with reference to the published genome 
sequence (<http://www.phytozome.net/mimulus>). 
Nested fragments of MgPDS (Fig. 1), MgCYCI, and 
MgCYC2 were amplified with primers PDS-F-BamHI 
(5’-ATAGGATCCTTGCW ATGCCAARY AARCC- 
WGG-3’) and PDS-R-Xhol (5’-AGTCTCGAGA- 
CACTIAGMAGKGRRCTTCTGC-3’), MgCYC1. 
VIGS.BamHIF (5’-ATAGGATCCGCAGCAAAGTAA- 
TAACCCTTCG-3’) and MgCYC1.VIGS.XhoIR_ (5/- 
AGTCTCGAGTCCTTCGAACCAGCTCCTTTAC-3"), 
and MgCYC2.VIGS.KpnIF (5’/-AGTGGTACCCGG- 
CATTGCTAGAAAGTTCTTC-3’) and MgCYC2. full. R 
to generate 5’ and 3’ restriction enzyme recognition 
sites. The resulting polymerase chain reaction (PCR) 
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Figure 1. Example VIGS constructs. —A. Structure of MgPDS showing the 312 bp fragment (black) cloned into TRV2, the 
5’ (gray) and 3’ (white) ends, and the position of the stop codon (star). Primer pairs used for amplification and sequencing are 
shown with arrows and dashed lines. —B. Map of the MgPDS-TRV2 construct used for VIGS. —C. Map of the Empty-TRV2 
construct used as an experimental control. LB, left border; RB, right border; Rz, self-cleaving ribozyme; NOSt, NOS terminator. 


fragments were digested with BamHI and Xhol (PDS 
and MgCYC1), or Kpnl and Xhol (MgCYC2) and 
subcloned individually or together (MgCYCI and 
MgCYC2) into the tobacco rattle virus 2 (TRV2) 


vector digested with the same enzymes (Fig. 1). 


ELECTROPORATION OF AGROBACTERIUM TUMEFACIENS WITH 
MODIFIED TRV2 


To optimize infection of Mimulus guttatus, MgPDS- 
TRV2, or Empty-TRV2, and TRV1 were used to 
electroporate two strains of Agrobacterium tumefaciens 
Smith & Townsen (GV3101 and EHA105) (Dinesh- 
Kumar et al., 2003; Hileman et al., 2005). Agro- 
bacterium tumefaciens was prepared according to a 
modified protocol (Hileman et al., 2005). Single 
colonies were used to inoculate 5 mL Luria Broth (LB) 
and were screened for the presence of the appropriate 
plasmids using primers pYLI56F/pYLI56R (TRV2 
constructs) and OYL195/OYL198 (TRV1) as previ- 
ously described (Hileman et al., 2005). PCR-positive 
cultures were used to inoculate 500 mL LB and were 
grown to an absorbance at 260 nm (Abs269) of seven 
to nine. Pelleted cultures were resuspended in 
infiltration media (Hileman et al., 2005) to a final 
Absv69 of two, incubated for four to 12 hr. at room 
temperature, and used to inoculate M. guttatus PR 
and IM767 seedlings. After optimization of VIGS, 
MgCYCI1-TRV2, MgCYC2-TRV2, and MgCYC1:Mg- 


CYC2-TRV2 were used to electroporate EHA105 as 
described. 


SEEDLING INOCULATION WITH VIGS CONSTRUCTS 


Mimulus guttatus seedlings were grown to the two 
to three leaf pair stage at 20°C-—22°C under long-day 
(16 hr. light/8 hr. dark) conditions on a 1:1 mix of 
vermiculite and soil. One hundred to 200 seedlings 
were inoculated with a 1:1 ratio of Agrobacterium 
tumefaciens GV3101 or EHA105 carrying TRV1 and 
MgPDS-TRV2, or Empty-TRV2 per treatment as 
previously described (Hileman et al., 2005). Batches 
of 50 M. guttatus PR and IM767 seedlings were 
subjected to four different inoculation treatments: 1 
min. vacuum infiltration, 3 min. vacuum infiltration, 
15 min. submersion without vacuum, and injection of 
leaves with a needleless syringe. Inoculated seedlings 
were replanted in soil and grown for an additional two 
(IM767) to four (PR) weeks at 20°C—22°C before leaf 
tissue was collected for RNA extraction. Differences 
in tissue harvest times were based on interpopulation 
variation in phenotype emergence. All surviving 
plants were examined for evidence of photobleaching 
as well as any obvious differences in flowering time 
and flower development relative to wild type. Once 
conditions had been optimized, VIGS was carried out 
as previously described in approximately 200 IM767 
seedlings using the MgCYCI-TRV2, MgCYC2- 
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TRV2, or MgCYC]:MgCYC2-TRV2 construct with 1 


min. vacuum infiltration. 


RNA EXTRACTION AND GENE EXPRESSION ANALYSIS 


Leaf tissue from surviving plants was harvested 
and pooled per five to 10 individual plants, and total 
RNA was extracted using TriReagent (Applied 
Biosystems, Foster City, California). For each RNA 
pool, 1 we was reverse transcribed using the iScript 
cDNA Synthesis Kit (Bio-Rad). cDNA was diluted 
1:10 and 2 uL was used as a PCR template to screen 
for the presence of TRV2 and TRV1 plasmids as 
previously described (Hileman et al., 2005). After 
identification of infected pools, total RNA was 
extracted from leaves of individual plants. cDNA 
was then synthesized after DNA digestion with 
TURBO DNA-free (Applied Biosystems), and each 
cDNA was screened for presence of TRV2 and 
TRV1. ACTIN was amplified for 28 cycles as an 
internal control as described in Prasad et al. (2001). 

To determine if a reduction in PDS levels was 
correlated with photobleaching, a ca. 150 bp region of 
PDS was amplified from cDNA of green leaf RNA 
harvested from Empty-TRV2-infected and uninfected 
individuals, and white leaf RNA harvested from 
MgPDS-TRV2-infected individuals. MgPDSrt-F (5/- 
CGAAGATGGCATTTTTGGAT-3’) and MgPDSrt-R 
(5’-ATGGCATCTCTTTCCATTCG-3’) were designed 
in Primer3 (Rozen & Skaletsky, 2000) to amplify a 
region downstream of the VIGS PDS insert region, 
and PCRs were run for 28, 30, and 32 cycles to 
confirm linearity. MgPDS expression levels were 
normalized against ACTIN in ImageJ (Abramoff et al., 
2004) with two to six biological replicates. 

Silencing of MgCYCI] and MgCYC2 was determined 
using quantitative (q) reverse transcriptase (RT)-PCR. 
The MgCYCI primers MgCYC1.rt.new.F (5'-CGAA- 
CAACGGGAAATTATGC-3') and MgCYC1.rt.new.R 
(5’-GCCGTAGTTGGAACTCGAAA-3'), and the 
MgCYC2 primers MgCY C2.rt.new.F (5’-CGGCCCGT- 
CTGCTAATAAT-3’) and MgCYC2.rt.new.R (5’- 
CCTATTGATGAACTTGGCTGGT-3’) were designed 
to amplify a 150 bp fragment downstream of the VIGS 
construct sequence in Primer3 (Rozen & Skaletsky, 
2000). Primer efficiency was determined using Fast 
SYBR Green Master Mix (Life Technologies) as 
previously described (Preston & Hileman, 2010). 
After correcting for transcriptional stability, cycle 
threshold values were normalized against the house- 
keeping gene EFlalpha as previously described 
(Scoville et al., 2011), and the mean was calculated 
for four technical replicates. Flowers positive for 


MeCYCI-TRV2, MgCYC2-TRV2, and MgCYC1:Mg- 
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CYC2-TRV2 were examined for organ number, shape, 
and identity, and compared to wild type flowers. 


RESULTS AND DISCUSSION 
PHOTOBLEACHING IN MGPDS-TRV2, VIGS TREATED PLANTS 


Silencing of PDS causes a characteristic photo- 
bleached phenotype that can be easily screened by 
the naked eye (Benedito et al., 2004; Hileman et al., 
2005). To compare the efficiency of VIGS for the two 
populations of Mimulus guttatus, we attempted to 
infect plants with two Agrobacterium tumefaciens 
strains, using four different treatments. Each 
MgPDS-TRV2-treated plant was then screened for 
evidence of photobleaching. For the PR population, 
MgPDS-TRV2-treated seedling survival ranged from 
68% to 84% using GV3101 and 3% to 29% using 
EHA105, with the highest survival rate for the 3 min. 
vacuum infiltration treatment for both A. tumefaciens 
strains (Table 1). Photobleached plants were identi- 
fied in all but the leaf injection treatment for both 
GV3101 and EHAI105. For the three successful 
treatments using GV3101 (1%-6%), the percentage 
of photobleached survivors was significantly higher 
(6%) after 3 min. vacuum infiltration (Table 1). This 
trend was also observed using EHA105. However, 
treatments using EHA105 were less efficient (0O%- 
3% infected) than GV3101 (Table 1). Furthermore, 
for all plants showing a phenotype, photobleaching 
was only observed on one or two leaves (Fig. 2A, B); 
no photobleaching was observed in reproductive 
tissues. 

Similar to the PR population, no photobleaching 
was observed for the leaf injection treatment in 
IM767. Furthermore, unlike the PR population, 
IM767 seedlings were unable to survive 3 min. 
vacuum infiltration regardless of the Agrobacterium 
strain used. For the other two treatments, survival 
ranged from 25% (1 min. vacuum; EHA105) to 40% 
(15 min. dip; both Agrobacterium strains). However, 
photobleached individuals were only found after 1 
min. vacuum infiltration. The efficiency of MgPDS- 
TRV2 VIGS using GV3101 (3%) was similar to the 
most efficient treatments in the PR population. By 
contrast, | min. vacuum infiltration using EHA105 
was significantly higher in the IM767 population 
(24%) versus the PR population (Table 1). This level 
of VIGS efficiency is similar to that found for PDS 
silencing in Papaver L. using EHA105 (Hileman et 
al., 2005) and Aquilegia L. using GV3101 (Gould & 
Kramer, 2007); VIGS has been successfully used in 
both these species to elucidate the function of flower 
developmental genes (Drea et al., 2007; Kramer et 


al., 2007; Yellina et al., 2010; Hands et al., 2011). 
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Table 1. Survival and infection rates for VIGS-treated Mimulus guttatus DC. seedlings. 
Point Reyes Iron Mountain 767 
Agrobacterium strain l min. 3 min. 3 min. 15 min. l min l min. 3 min. 15 min. 
GV3101 
Construct PDS E PDS PDS E PDS PDS PDS 
N treated 150 100 100 200 100 100 100 100 
Survival (%) 68 60 84 67 40 34 0 40 
Survivors infected (%) 2 0 6 1 0 3 N/A 0 
Phenotype! (%) 2 0 6 l 0 3 N/A 0 
EHA105 
Construct PDS E PDS PDS E PDS PDS PDS 
N treated 150 100 100 150 100 150 100 100 
Survival (%) 3 24 29 24 iA 25 0 40 
Survivors infected (%) 0 0 3 0 14 24 N/A 0 
Phenotype! (%) 0 0 3 0 0 24 N/A 0 


Abbreviations: PDS, MgPDS-TRV2; E, Empty-TRV2. 


! Percentage of survivors. 


Furthermore, ca. 50% of the IM767 plants showing a 
phenotype were almost completely photobleached 
(Fig. 2D, E). This level of photobleaching caused 
severe retardation (Fig. 2D, E), followed by death 
prior to development of the reproductive organs. 
Consequently, no photobleaching was observed in 
reproductive structures. 

To confirm that photobleaching in both Mimulus 
guttatus populations was positively correlated with 
MgPDS-TRV2 and TRV1 infection, and was not a 
result of infection per se, tissues from both MgPDS- 
TRV2- and Empty-TRV2-treated plants were PCR 
screened for the presence of TRV2 and TRV1. For 
both the PR and IM767 populations, all photo- 
bleached leaves were infected with both MgPDS- 
TRV2 and TRVI1 (Fig. 3A); no green leaves were 
positive for either vector. Control plants positive for 
both Empty-TRV2 and TRV1 were only found in the 
IM767 population using EHA105 (14%), consistent 
with infection rates for MgPDS-TRV2 (Table 1; Fig. 
3A). None of these infected control plants showed 
photobleaching in vegetative or reproductive tissues, 
and growth was normal throughout development. 
Thus, photobleaching was correlated with MgPDS- 
TRV2 infection, rather than viral infection per se. 
Furthermore, all Empty-TRV2-treated plants that 
screened positive for the virus in leaves also screened 
positive in flowers, suggesting that infection remains 
into the reproductive stage of development. 


PDS TRANSCRIPT LEVELS IN PHOTOBLEACHED MIMULUS 
GUTTATUS PLANTS 


To verify that photobleaching in MgPDS-TRV2- 
infected plants correlated with reduced endogenous 


MgPDS transcript levels, we performed RT-PCR on 


both PR and IM767 plants infected with MgPDS- 
TRV2 in GV3101 and EHA105, IM767 plants 
infected with Empty-TRV2 in EHA105, and plants 
from the same treatments that were not infected. For 
each infected plant, endogenous MgPDS expression 
levels were calculated relative to uninfected plants 
from the same treatment (set to an expression level of 
1) after standardization against ACTIN (Fig. 3B). RT- 
PCR analyses confirmed silencing of endogenous 
MgPDS in MgPDS-TRV2 plants relative to both 
uninfected and Empty-TRV2 control plants (Fig. 3B). 
For both GV3101 and EHA105 treatments, average 
silencing was two-fold with a high level of variation 
between individuals, ranging from no significant 
silencing to 10-fold silencing (Fig. 3B). Because 
whole leaves were harvested for expression analyses, 
inter-individual variation in endogenous MgPDS 
expression is most likely due to differences in the 
number of leaf cells infected with MgPDS-TRV2. 
The amount of photobleaching per MgPDS-TRV2- 
infected individual ranged from small sectors on 
individual leaves (GV3101 in PR) to complete 
photobleaching of one or two leaves (EHA105 in 
PR and IM767) to almost complete photobleaching of 
the whole seedling (EHA105 and GV3101 in IM767) 
(Fig. 2A, B, D, E). This level of VIGS photobleaching 
variation is consistent with reports from other species 


(Hileman et al., 2005; Liu & Page, 2008). 


TRANSCRIPT SILENCING OF PUTATIVE FLORAL CYC SYMMETRY 
GENES 


The fact that Mimulus guttatus IM767 infection 
persists into flowers (previous section) suggests 
efficacy of VIGS in reproductive tissues. However, 
in order to verify that infection of reproductive tissues 


306 Annals of the 
Missouri Botanical Garden 


Figure 2. Variation in leaf photobleaching and flower morphology for MgPDS-TRV2- and MgCYC-TRV2-infected plants. A, 
B. Infected plants of Mimulus guttatus PR treated with MgPDS-TRV2 in GV3101 show variation in the extent of photobleaching 
[arrow in (B)| in individual leaves. —-C. Uninfected M. guttatus PR plants treated with MgPDS-TRV2 show no sign of 
photobleaching. —D. Almost completely photobleached MgPDS-TRV2 in EHA105-infected M. guttatus IM767 seedling. —E. 
Initial stages of systemic photobleaching for an MgPDS-TRV2 in GV3101-infected M. guttatus IM767 seedling. —F. 
Uninfected M. guttatus IM767 plants treated with MgPDS-TRV2 show no sign of photobleaching. —G. Empty-TRV2 positive 
flower showing wild type morphology. —H. MgCYC1-TRV2-infected flowers (right) have small dorsal petal and wavy petal 
margins compared to uninfected flowers (left). —I. Infected flower in H. Infected flower prior to anthesis has shortened dorsal 
petals, resulting in protrusion of the stigma. —J. Another example of an MgCYC1-TRV2-infected flower showing wavy margins 
in all petals. K-M. MgCYCI:MgCYC2-TRV2-infected flowers have an abnormal number of lateral and/or ventral 
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also causes gene silencing and has a correlated effect 
on phenotype, we silenced two homologs of the 
Antirrhinum majus clade Il TCP family flower 
symmetry genes CYCLOIDEA (CYC) and DICHOTO- 
MA (DICH) (Luo et al., 1996, 1999). In A. majus, 
mutations in CYC result in partial loss of dorsal petal 
and staminode identity, resulting in plants with five 
fully developed stamens and partially radialized 
petals (Luo et al., 1996, 1999). By contrast, the dich 
mutant has no effect on stamen number but causes 
the loss of internal asymmetry within the dorsal 
petals; loss-of-function mutations in both CYC and 
DICH result in fully radialized flowers that have six 
ventral petals and six ventral stamens (Luo et al., 
1996, 1999). 

Of the IM767 seedlings vacuum infiltrated for 1 
min. with the MgCYC gene-silencing constructs, 
three (6% of survivors), five (28% of survivors), and 
six (15% of survivors) individuals had flowers that 
were positive for MgCYCI-TRV2, MgCYC2-TRV2, 
and MgCYC1:MgCYC2-TRV2, respectively. Q-RT- 
PCR analyses of gene expression revealed that 
MgCYCI expression was significantly reduced in 
MgCYCI-TRV2 and MgCYCI:MgCYC2-TRV2- 
infected versus uninfected flowers collected for 
RNA post-anthesis (Fig. 3C). Similarly, MgCYC2 
expression was significantly reduced in MgCYC2- 
TRV2 and MgCYC1:MgCYC2-TRV2-infected versus 
uninfected post-anthesis flowers (Fig. 3C). This 
demonstrates that gene silencing was target specific. 
In the case of MgC YC2-silenced flowers, no abnormal 
phenotype was detected relative to uninfected plants 
in terms of organ number, size, shape, or identity. By 
contrast, infection and partial silencing of MgCYC1, 
and MgCYC1 and MgCYC2 together was correlated 
with abnormal petal and/or stamen whorl phenotypes 
but did not affect sepal morphology (Fig. 2H—O). 

All nine MgCYCI-TRV2-infected flowers from 
three plants had petals with wavy margins compared 
to the smooth margins of wild type and Empty-TRV2- 
infected flowers (Fig. 2G—J). In each case, this growth 
defect was found in dorsal, lateral, and ventral petals 
(Fig. 2H-J). In addition to wavy petal margins, the 
two dorsal petals of infected flowers appeared 
shortened relative to wild type dorsal petals and 
infected lateroventral petals (Fig. 2H—J). Quantitative 
analyses of flower size revealed that MgCYC1-TRV2- 
infected flowers were shorter and narrower than 
uninfected flowers (Fig. 3D, E). Furthermore, the 
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average ratio between dorsal and ventral petal lobe 
length was smaller in MgCYC1-TRV2-infected 
relative to uninfected flowers (Fig. 3F). This suggests 
that the reduction in growth was more extreme for 
dorsal petals than for ventral petals. Indeed, in some 
cases, reduction of dorsal petal size was evident at the 
bud stage (Fig. 2I); in wild type flower buds, the 
dorsal petals enclose the reproductive structures, 
whereas in two MgCYC1-TRV2-infected flowers, the 
anthers and stigma protruded from the flower bud. 

In striking contrast to MgCYCI- and MgCYC2- 
silenced flowers, all MgCYC1:MgCYC2-TRV2-infect- 
ed flowers had extra petals and/or stamens relative to 
wild type and Empty-TRV2-infected flowers (Figs. 
2L—O, 3G). In wild type Mimulus guttatus plants, five 
stamen primordia give rise to two lateral and two 
ventral stamens and a single dorsal staminode that 
arrests growth in early development (pers. obs.). By 
contrast, in MgCYC1:MgCYC2-TRV2-infected flow- 
ers, one or two stamens developed in the dorsal 
region, resulting in a total of five to six pollen- 
producing stamens at anthesis (Figs. 2M, 3G). This 
gene-silenced phenotype is similar to cyc and 
cyc:dich mutants of Antirrhinum majus, which 
develop five and six stamens, respectively, compared 
to four stamens and one dorsal staminode of wild type 
flowers (Luo et al., 1996, 1999). In addition to 
stamens, MgCYC1:MgCYC2-TRV2-infected flowers 
showed organ number and identity defects in the 
petal whorl. Infected flowers ranged from having one 
less lateral petal to having both extra lateral and 
ventral petals (Figs. 2K—O, 3G). In the most extreme 
case, flowers were almost fully radialized due to the 
development of three sets of two lateral and one 
ventral petal. The flowers of this plant also had six 
stamens and showed a novel phenotype consisting of 
two instead of one gynoecium (Fig. 2N, O). 


DIVERGENCE OF CYC GENE FUNCTION IN ASTERIDS 


MgCYC1I and MgCYC2 are the closest homologs of 
CYC and DICH from Antirrhinum majus, falling into 
the ECE CYC2 clade of TCP genes (Howarth & 
Donoghue, 2006). However, the duplication events 
that gave rise to the ECE CYC2 clade paralogs of 
Mimulus guttatus and A. majus occurred indepen- 
dently. Thus, comparison of these genes in asterids 
can inform our understanding of functional diversi- 
fication after gene duplication. The VIGS data 


petals (asterisks). —N. Comparison of an extreme MgCYC1:MgCYC2-TRV2-silenced floral phenotype (left) with that of an 
uninfected flower (right). —O. Dissected flower from (N) showing an increase in stamen, petal, and gynoecium number. In 
addition to a doubling of the lateroventral petal unit in the ventral domain, the dorsal petals are partially transformed into a 


dorsal lateroventral-like unit. D, dorsal; L, lateral; V, ventral. 
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Figure 3. VIGS-mediated silencing of PDS and CYC genes. —A. RT-PCR showing the presence of TRV vectors in two 
representative IM767 individuals infected with MgPDS-TRV2 (PDS; larger TRV2 band) and Empty-TRV2 (E; smaller TRV2 
band). These screening methods were used to determine the number of infected individuals described in Table 1. Endogenous 
levels of MgPDS are reduced in individuals infected with MgPDS-TRV2 relative to Empty-TRV2-infected and uninfected 
individuals. ACTIN was amplified as an internal control. —B. Boxplots showing variation in ACTIN-standardized endogenous 
MgPDS expression levels relative to the average of two to four uninfected individuals (dotted line) for each population and 
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presented in this paper support both similar and 
divergent functions for MgCYCI and MgCYC2 in 
flower development and suggest differential sub- and 
neo-functionalization of paralogous genes after inde- 
pendent duplications in Plantaginaceae (A. majus) 
and Phrymaceae (M. guttatus). 

In the case of MgCYC1, silencing has an effect on 
dorsal petal size, consistent with the known function 
of TCP genes on dorsal petal development and cell 
proliferation (Luo et al., 1996, 1999; Feng et al., 
2006; Busch & Zachgo, 2007; Broholm et al., 2008; 
Preston & Hileman, 2009; Kieffer et al., 2011; Bai et 
al., 2012; Braun et al., 2012). However, within the 
ECE CYC2 clade, the specific phenotype observed in 
MgCYCl-silenced flowers is unique to Mimulus 
guttatus. In Antirrhinum majus, mutations in either 
CYC or DICH have negative effects on dorsal petal 
growth, but mutations in CYC have positive effects on 
stamen growth. Furthermore, lateral petals of double 
cyc:dich mutants are transformed into ventral petals, 
and both petal and stamen number is increased (Luo 
et al., 1996, 1999). A role in dorsal identity is also 
suggested for species in the Commelinaceae, Poa- 
ceae, Malpighiaceae, Caprifoliaceae, Gesneriaceae, 
and Papaveraceae based on gene expression data 
(Damerval et al., 2007; Song et al., 2009; Yuan et al., 
2009; Zhang et al., 2010; Howarth et al., 2011; 
Preston & Hileman, 2012). 

Similar to Antirrhinum majus cyc mutants, silenc- 
ing of MgCYCI negatively affects dorsal petal growth. 
However, like A. majus dich mutants, Mimulus 
guttatus dorsal petals are not transformed into lateral 
or ventral petals, and there is no effect on floral organ 
number. MgCYCI silencing also causes growth 
defects at the margins of dorsal, lateral, and ventral 
petals; this phenotype is so far novel to loss-of- 
function ECE CYC2 clade genes. Interestingly, a 
similar phenotype was recently observed when 
several distantly related CINCINATTA-like (CIN- 
like) TCP genes were simultaneously silenced 
(Koyama et al., 2011). Given the lack of this 
phenotype in Plantaginaceae and Asteraceae CYC- 
like mutant backgrounds, it is most parsimonious to 
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infer the recent acquisition of petal margin growth 
function in the MgCYCI lineage. However, it is 
unclear why this phenotype was not observed in 
MgCYC1:MgCYC2-TRV2-infected plants. Our work- 
ing hypothesis is that the loss of competition between 
MgCYCl1 and MgCYC2 for protein partners may have 
resulted in a gain of function phenotype; this remains 
to be tested. 

Although silencing of MgCYC2 in a few Mimulus 
guttatus flowers resulted in no abnormal phenotypes, 
the observation that double MgCYCI:MgCYC2- 
TRV2-infected flowers on average had more petals 
and stamens than single gene-silenced flowers 
suggests that MgCYCI and MgCYC2 have a 
redundant role in specifying dorsal identity and 
stabilizing organ number. A similar role has been 
assigned to CYC and DICH and is suggested for other 
asterid CYC-like genes (Luo et al., 1996, 1999; 
Cubas et al., 1999; Gao et al., 2008; Song et al., 
2009; Preston et al., 2009, 2011; Fambrini et al., 
2011). However, the fact that MgCYC1:MgCYC2- 
TRV2-infected flowers tend to duplicate the later- 
oventral petals as a single unit tentatively suggests 
differential integration of the petal module in 
asterids. This hypothesis warrants further investiga- 
tion in the future. 


CONCLUSION 


By applying different treatments to two ecologically 
and morphologically distinct populations, we have 
optimized VIGS as a reverse genetics tool to assay 
protein function in Mimulus guttatus. Specifically, 
Agrobacterium tumefaciens EHA105-mediated treat- 
ments of IM767 seedlings at the two to three leaf pair 
stage, via l min. vacuum infiltration, proved highly 
efficient. This treatment of IM767 seedlings yielded 
ca. 14%-24% infected survivors, the majority of 
which carried the virus in all tissues throughout 
development (Table 1). In the case of infection with 
MgPDS-TRV2, systemic infection caused severe 
photobleaching, resulting in death prior to flowering. 
Thus, it was not possible to determine the ability of 
the virus to infect reproductive tissues. By contrast, 


treatment. Endogenous MgPDS levels are reduced in MgPDS-TRV2-infected individuals, but not Empty-TRV2-infected 
controls, relative to uninfected individuals using both EHA105 and GV3101 Agrobacterium strains. Data for MgPDS-TRV2- 
infected individuals include both the PR and IM767 populations; Empty-TRV2-infected individuals are from the IM767 
population only. Solid horizontal line, median; box, interquartile range; whiskers, extending to the further point that is within 
1.5X the interquartile range from the box. —C. Quantitative-RT-PCR expression data for MgCYC1 (black bars) and MgCYC2 
(gray bars) in uninfected and MgCYCI:MgCYC2-TRV2-infected flowers post-anthesis. Data are averages for three biological 
replicates. Error bars denote standard deviations. —D. Effect of CYC] silencing on flower width. —E. Effect of CYCT silencing 
on flower height. —F. Effect of CYC silencing on the ratio of dorsal to ventral (DV) petal lobe length. —G. Effect of 
simultaneous CYC] and CYC2 silencing on organ number. Lateral petal (LP), ventral petal (VP), stamen (St), and gynoecium 
(Gy), but not dorsal petal (DP), numbers are all increased relative to uninfected control plants (dotted line). Error bars denote 


standard deviations. 
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screening for both TRV2 and TRV1 in individuals 
infected with the Empty-TRV2 control vector re- 
vealed 100% retention of the virus in flowers, without 
any effect on flowering time, vegetative and inflores- 
cence development, or seed set. The efficacy of VIGS 
in flowers was verified by the silencing of two CYC- 
like homologs, resulting in wavy petal margins, 
reduced dorsal petal growth, extra lateroventral petals 
and stamens, and variable loss of dorsal identity. 
Together, these data demonstrate the potential of 
VIGS as a fast and reliable method to assay gene 
function in M. guttatus, advancing this species 
further as a tool for addressing fundamental ecolog- 
ical and evolutionary questions. 
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